Study of the ao(980) meson via the radiative decay 
(f) r]7T^j with the KLOE detector 
KLOE Collaboration 



F. Ambrosino° ^ A. Antonelli^, M. Antonelli^, F. Archillij '^, C. Bacci^ ™, 
P. Beltrame'', G. Bencivenni'', S. Bertolucci'', C. Bini'^ ', C. Bloise^, 
S. Bocchetta' ™, F. Bossi'', P. Branchini™, P. Campana'', G. Capon*^, 
T. Capussela^, F. Ceradini''", S. ChP, G. Chiefari^f, P. Ciambrone^, 
F. Crucianelli^, E. De Lucia'', A. De Santis^ ', P. De Simone'^, G. De Zorzi'^ ', 
A. Denig^, A. Di Domenico'^ ', C. Di Donato^, B. Di Micco' ™, A. Doria^, 
M. Dreucci^, G. Felici^, A. Ferrari'*, M. L. Ferrer'*, S. Fiore^-', C. Forti'*, 
P. Franzini'* ', C. Gatti'*, P. Gauzzi*'^ ', S. Giovannella'*, E. Gorini'^ 'i, 
E. Graziani", W. Kluge'^, V. KulikovP, F. Lacava^ ', G. Lanfranchi'*, 
J. Lee-Franzini'*-", D. Leone**, M. Martini'*s, P. Massarotti^' f, W. Met'*, 
S. Meola'^-^ S. Miscetti'*, M. Moulson'*, S. Miiller^, F. Murtas'*, 
M. Napolitano^-^ F. Nguyen''™, M. Palutan'*, E. Pasqualucci', A. Passeri™, 
V. Patera^^s, F. Perfetto'=f, M. Primavera'^, P. Santangelo'*, G. Saracino°^ 

B. Sciascia'*, A. Sciubba'*s, A. Sibidanov'*, T. Spadaro'*, M. Testa'* ', 
L. Tortora***, P. Valente', G. Venanzoni'*, R.Versaci'*^, R. Volpe'* ', G. Xu'*" 

"■Laboratori Nazionali di Frascati deU'INFN, Via E.Fermi 4O, 1-00044 Frascati, Italy. 
Institut fiir Experimentelle Kernphysik, Universitdt Karlsruhe D-76128 Karlsruhe, 

Germany. 

Dipartimento di Fisica dell'Universitd del Salento, Via Arnesano, 1-73100 Lecce, Italy. 
'^INFN Sezione di Lecce, Via Arnesano, 1-73100 Lecce, Italy. 
" Dipartimento di Scienze Fisiche deW Universita "Federico 11", Via Cintia, 1-80126 Napoli, 

Italy 

fiNFN Sezione di Napoli, Via Cintia, 1-80126 Napoli, Italy 
^Dipartimento di Energetica, Sapienza Universita di Roma, P.le A.Moro, 2 1-00185 Roma, 

Italy. 

'^Dipartimento di Fisica, Sapienza Universita di Roma, P.le A.Moro, 2 1-00185 Roma, Italy. 
'INFN Sezione di Roma, P.le A.Moro, 2 1-00185 Roma, Italy. 
^Dipartimento di Fisica dell' Universita di Roma "Tor Vergata" , Via della Ricerca 
Scientifica, 1 1-00133 Roma, Italy. 
^INFN Sezione di Roma Tor Vergata, Via della Ricerca Scientifica, 1 1-00133 Roma, Italy. 
^Dipartimento di Fisica dell'Universitd di Roma "Roma Tre", Via della Vasca Navale, 84 

I-OOI46 Roma, Italy. 
"^INFN Sezione di Roma Tre, Via della Vasca Navale, 84 I-OOI46 Roma, Italy. 
"Physics Department, State University of New York at Stony Brook, Stony Brook, NY 

11794-3840 USA. 

"Institute of High Energy Physics of Academica Sinica, PO Box 918 Beijing 100049, PR 

China. 

f Institute for Theoretical and Experimental Physics, B. Cheremushkinskaya ul. 25, 
RU-117218 Moscow, Russia. 



* Corresponding author, Dipartimento di Fisica, Sapienza Universita di Roma, P.le A.Moro, 
2 1-00185 Roma, Italy; Tel. +390649914266, Fax +39064957697 
Email address: paolo.gauzziaromal.infn.it (P. Gauzzi) 



Preprint submitted to Physics Letters B 



April 21, 2009 



Abstract 



We have studied the oq (980)7 process with the KLOE detector at the 
Frascati (/)— factory DA$NE by detecting the (j) rjir'^j decays in the final 
states with rj — > 77 and 77 — s- tt+tt^tt'^. We have measured the branching 
ratios for both final states: Br{(t> 7771%) = (7.01 ± 0.10 ± 0.20) x 10"^ and 
(7.12 ± 0.13 ± 0.22) X 10"^ respectively. We have also extracted the ao(980) 
mass and its couplings to rjir'^, K^K~, and to the (j) meson from the fit of the 
rjTT^ invariant mass distributions using different phenomenological models. 

Key words: e~ collisions, Scalar mesons, Rare decays 
PACS: 12.39.Mk, 13.20.Jf, 13.66.Bc, 14.40.Cs 



1. Introduction 

The problem of the internal structure of the scalar mesons with mass below 
1 GeV is still open[l|. It is controversial whether they are qq mesons fij, qqqq 
statesQ, bound states of a KK pair[3| or a mixing of these configurations. 
An important part of the program of the KLOE experiment, carried out at the 
Frascati ^-factory DA<I>NE, has been dedicated to the study of the radiative 
decays ^(1020) P1P2I (A, 2 = pseudoscalar mesons). These decays are 
dominated by the exchange of a scalar meson S in the intermediate state {(j) — » 
5*7, and S — *■ P1P2), and both their branching ratios and the P1P2 invariant 
mass shapes depend on the scalar structure. 

The (p 777r°7 decay has been already used by KLOE and by other experiments 
to study the neutral component of the isotriplet ao(98O)0, Q. This process is 
well suited to study the (j) — > oq (980)7 dynamics, since it is dominated by the 
scalar production, with small vector background, contrary to tt'^tt^j and 7r+7r~7 
cases, where a large irreducible background interferes with the /o(980) signalQ- 
In this paper the result of the analysis of the (j) rjir^j decay, performed on 
a sample with 20 times larger statistics than the previously published paper 
is presented. The final states corresponding to 77 — *■ 77 and ij —> Tr+Tr^Tr" have 
been selected. The ijir^ invariant mass distributions have been fit to two models 
of parametrization of the oq (980)7 decay, to extract the relevant ao(980) 
parameters (mass and couplings). 

2. DA*NE and KLOE 

The Frascati 0-factory DA<i>NE is an e+e^ collider operating at a center 
of mass energy y/s = ~ 1020 MeV. The beams collide at an angle of (tt - 
0.025) rad, thus producing (j) mesons with small momentum [p^ ~ 13 MeV) in 
the horizontal plane. The KLOE detector^S,] consists of two main subdetectors: 
a large volume drift chamber (DC) and a fine sampling lead-scintillating fibers 
electromagnetic calorimeter (EMC). The whole apparatus is inserted in a 0.52 T 
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axial magnetic field, produced by a superconducting coil. The DC is 3.3 m long, 
with inner and outer radii of 25 and 200 cm respectively. It contains 12 582 drift 
cells arranged in 58 stereo layers uniformly distributed in the sensitive volume 
and it is filled with a gas mixture of 90% helium and 10% isobutane. Its spatial 
resolution is 200 /zm and the tracks coming from the beam interaction point 
(IP) are reconstructed with (T{p±)/p± < 0.4%. The position resolution for two 
track vertices is about 3 mm. 

The DC is surrounded by the EMC, that covers 98% of the sohd angle, and is 
divided into a barrel, made of 24 trapezoidal modules about 4 m long, with the 
fibres running parallel to the barrel axis, and two endcaps of 32 module each, 
with fibers aligned vertically. The read-out granularity is ^ 4.4 x 4.4 cm^, for 
a total of 2440 cells, read at both ends by photomultiplicrs. The coordinate 
of a particle along the fiber direction is reconstructed from the difference of 
the arrival time of the signals at the two ends of the cell. Cells close in time 
and space are grouped together into clusters. The cluster energy is the sum of 
the cell energies, while the cluster time and position are energy weighed aver- 
ages. The energy and time resolutions for photons are aE 

and at — 57 ps/-\/i?(GeV) 100 ps respectively. Cluster positions are mea- 
sured with resolutions of 1.3 cm in the coordinates transverse to the fibers, and 
1.2 cm/ y/E{GeV) in the longitudinal coordinate. The detection efficiency for 
photons of £^ ~ 20 MeV is greater than 80% and reaches almost 100% at i? > 80 
MeV. 

The KLOE trigger is based on the detection of two energy deposits in the EMC, 
with E > 50 MeV in the barrel and E > 150 MeV in the endcaps. 

3. Event selection 

The results are based on the data collected during the 2001-02 run, at 
y/s ~ M^. Of the two selected decay chains, the fully neutral one is char- 
acterized by high statistics and large background, while the charged one has 
small background but lower statistics. These two decay chains have been se- 
lected with different criteria and slightly different data samples have been used: 
414 pb^^ for the fully neutral and 383 pb^^ for the charged decay. Monte Carlo 
(MC) samples of signal and of backgroundprocesses have been produced with 
the simulation program of the experiment [9|. They have been generated on a 
run-by-run basis, simulating the machine operating conditions and background 
levels as measured in the data. Three MC samples, generated with different 
luminosity scale factors (LSF = Lmc / Ldata), have been used: 

1. the rad sample contains all the radiative (/)— decays plus the non resonant 
process e^e~ — ^ cjtt*', with LSF=5; 

2. the kk sample contains </> K'^K^ with all subsequent kaon decays gen- 
erated with LSF=1; 

3. the all sample contains all the decays with LSF=l/5; it is used to find 
possible backgrounds not included in the two main samples. 
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The shape of the t^tt" invariant mass distribution has been simulated according 
to the spectrum obtained from the previously published analysis 

3.1. (/)—!■ ?77r°7 with rj ^ 

This final state is characterized by five prompt photons originating from the 
IP. A prompt photon is defined as an EMC cluster not associated to any charged 
track in the DC and satisfying the condition \t~r/c\ < mm[5at{E) , 2 ns], where 
t is the photon flight time, r is the corresponding path length, and c is the speed 
of light. Events with exactly five prompt clusters, with E > 3 MeV and polar 
angle i? > 21° with respect to the beam line, are selected. 

The main background originates from the other five photon final states, (j) ~* 
/o (980)7 7r'^7r'^7 and e~^e~ — > ujtt'^ — 7r''7r''7, and from the seven photon 
process, (p ^ rjj with rj —^ Stt*^, which can mimic five photon events due to either 
photon loss or cluster merging. Also the three photon final states, (j) vi''^^)! 
with 7y(7r^) 77, give a small contribution to the selected sample, when fake 
clusters are produced either by accidental coincidence with machine background 
or by cluster splittings. Other background processes are negligible. 
The following analysis steps are then applied to the selected events. 

1. First kinematic fit which imposes the total 4- momentum conservation and 
the speed of light for each photon, with 9 degrees of freedom. Events 
with X/iti > 27 are rejected. A cut at 980 MeV on the total energy of 
the three most energetic photons is also applied to reject residual three 
photon events (processes 4 and 5 of Table [T]) . 

2. Search for the best photon pairing to 77's and tt^'s, by choosing the combi- 
nation that minimizes the x^-like variable {i,j,k,l — l,...,5 are the photon 
indices): 

2 ^ {M,,-Mp,f {Mki-Mp,f 

"Kpair 2 2 

for both P1P2 = r]TT° (signal) or 7r°7r° (background) hypotheses, ctm^o 
and ctm are the width of the tt*^ and rj peaks after the first kinematic fit 
(aM^^ = 6 MeV and = 9 MeV). 

3. Second kinematic fit with the two additional constraints of the masses of 
the intermediate particles. The number of degrees of freedom is 11. 

Background from process 1 and 3 of Table [T] dominates the tail of the distribu- 
tion of the Xfit2 of the second kinematic fit, as shown in Fig[Tl and it can be 
reduced by cutting at X%t2 ^ ^4. By using the photon pairing in the background 
hypothesis, 7r°7r°7, the Dahtz plot of FigUis obtained: the /o7 background pop- 
ulates the lower right corner, while the two straight bands are the contribution 
of ujtt'^. The oq signal is contained in the region between these bands. The utt^ 
background is strongly reduced by cutting out the two bands shown in FigU] 
Assuming the background hypotesis cjtt^, the angle 0* between the non as- 
sociated photon and the uj flight direction can be defined. The regions at 
large | cos 6**1 (Figl^left) are dominated by uj7t° and /o7 backgrounds. The 
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Figure 1: Right: of the second kinematic fit; the applied cut at Xjit2 ~ '^^ shown. 
Left: Dahtz plot of data in the background hypothesis (7r''7r''7). 



cut I COS 0*1 < 0.8 is then applied. Another effective cut to reduce the /o7 back- 
ground is 623 > 42° (FigHlright), where 623 is the angle between the second and 
third photons ordered by decreasing energy. 

After these cuts the overall selection efficiency, evaluated by MC, is almost in- 
dependent from the rjir'^ invariant mass and its average value is 38.5%. The final 




Figure 2: Left: cosS* distribution (see text for explanation); right: angle between the second 
and third photons ordered by decreasing energy (vertical lines represent the applied cuts). 



sample consists of 29 601 events and the expected S/B ratio is about 1.0 (see 
Table [1]). The residual background is irreducible and has to be evaluated and 
subtracted. A reweighing procedure has been adopted: 

1. for each specific background process a data sample with a small signal 
content (below few percent) has been selected; 

2. a fit has been performed on selected kinematical distributions, using the 
corresponding MC shapes to determine the weight to be assigned to that 
specific background; the weight is defined as the ratio of the number of 
events found by the fit to the number of expected events from MC. 
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In the last two columns of Table [T] the applied weights and the numbers of 
background events in the final sample are listed. The uncertainties are the com- 



Table 1: Background processes for (j) rjir'-^-y, with rj — + 77. (S/B)i is the signal to background 
ratio after the preselection, (S/B)2 the same ratio at the end of the whole analysis chain. The 
reweighing factors, w, are also listed. Last column reports the final background estimate. 





Process 


(S/B)i 


(S/B)2 


w 


Background events 


1 


/o7 ^ 7r'^7r"7 


0.40 


4.4 


1.2 


5062 ± 60 


2 




0.14 


3.1 


0.96 


3825 ± 37 


3 


4> ~^ VI with 77 — > Stt" 


0.10 


2.8 


1.1 


7248 ± 78 


4 


(f) — !• r/7 with 77 — !■ 77 


1.6 


200 


2.5 


197 ± 11 


5 


(p — !■ 7r*'7 


10 










Total background 


0.05 


1.0 




16 332 ± 86 



bination of MC statistics and of the systematics on the applied weights. The 
correlations have also been taken into account. After the background subtrac- 
tion the number of signal event is 13 269 ± 192. In Figl3]the ?y7r° invariant mass 
distribution of the final sample is shown together with the background contri- 
butions. The invariant mass resolution is about 4 MeV, with non-gaussian tails 
mainly due to wrong photon combinations. In the same figure, the distribution 
of the polar angle dree of the recoil photon is plotted after the background sub- 
traction: good agreement with the expected 1 -I- cos^ 9rec behaviour is shown. 




Figure 3: Left: jyTr" invariant mass distribution of the neutral channel. Right: Distribution 
of the cosine of the polar angle of the recoil photon after background subtraction (dots), 
compared with the MC expectation (solid line). 



3.2. (/)—*■ 7?7r°7 with rj — > tt~^7t tt*^ 

With respect to the fully neutral one, this decay provides a lower statistics 
since the branching ratio of 77 ^ tt+tt^tt'' is smaller than for 77 77. Moreover 
a lower acceptance is expected due to the larger number of particles to be 
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detected. However in this case there is a smaller background contamination, 
since no other final state with two tracks and five photons has a significant 
branching ratio from the (p. The main sources of background are due to final 
states with two tracks and either four or six photons. In order of importance 
there are: e+e~ — > wtt^ with w Tr+vr^Tr" and a fake cluster; (f) KsKl with 
Ks — > iT'^iT~ and prompt — > 377^^ with one photon lost; — > KsKl with 
Ks TT^TT^ and prompt Kl tt+tt^tt'^ or niiy with either one photon lost or 
one fake cluster; (j) —> rj^y with 77 — + tt^tt^tt'^ plus two fake clusters. 
The signal preselection requires the detection of two charged tracks and of five 
photons. The following requirements are then applied: 

1. a vertex with two opposite sign tracks in a cylinder, around the IP, of 5 
cm radius and 11 cm length; 

2. five prompt photons with E >10 MeV; 

3. total energy in the range 900 < Etot < 1160 MeV and total momentum 
\P7ot\ < llOMeV/c; 

4. the scalar sum of the momenta of the two pions Pj; = |pi | + |P2 1 , outside the 
range 418 < Ps < 430 MeV/c, which identifies events with Ks tt+tt^. 

Events surviving this preselection go to the kinematic fit stage, similar to that 
of the neutral channel. 

1. A kinematic fit with 9 degrees of freedom is performed by imposing only 
the total 4-momentum conservation and speed of light for the photons; 
events with xjui < ^'^ retained. 

2. Photons are combined to build tt^'s and rj's. There are 15 possibilities to 
get two 7r°'s out of five photons. For each of them there are two choices 
in the association of one 7r° to the tt+tt^ pair. For each of these 30 
combinations Xpair computed according to {i,j,k,l — 1,...,5 are the 
photon indices): 



2 _ (My - A4o)2 {Mm - M^of (M^+^-^o - 



Xpair 2 ' 2 ' 2 

Events with at least one combination with Xpair < 10 are retained. 

3. The second kinematic fit is performed on all the combinations selected 
by the previous step adding the three mass constraints, for a total of 12 
degrees of freedom. The combination with the lowest X/it2 chosen. Only 
events with X/it2 < ^0 are retained. 

4. Finally, events with the recoil photon energy below 20 MeV are discarded 
to remove events with a spurious low energy photon. 

The final sample consists of 4181 events. The overall selection efficiency for 
the signal, evaluated by MC, is 19.4%, almost independent from the rjir^ invari- 
ant mass, decreasing only at very high invariant mass values. Fig|4] shows the 
data-MC agreement for the distributions of the first and second kinematic 
fits. The MC distributions include signal and background events. The mass 
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Figure 4: distributions for the first (left) and second (right) kinematic fit. The selected 
data sample (points) is compared to the MC expectation (dark grey histograms) given by the 
weighed sum of the signal and the estimated background (light grey histograms). 



resolution is about 4 MeV for all mass values, with non gaussian tails, mainly 
due to events with a wrong photon combination. 

The residual background is evaluated by applying the selection procedure on 
MC samples and by checking the absolute normalization on background enriched 
data control samples. In order to properly normalize the observed numbers of 
events, data and MC samples after the preselection but before the kinematic 
fit have been used. At this level the expected contribution of the signal does 
not exceed 2 3%. Four variables have been chosen to compare data and MC 
samples: Etot, \Ptot\, M-y-y and M^7r77 where M-y^ is the invariant mass of any 
pair of photons (10 combinations per event) and M^^t^^^ is the invariant mass of 
the two pions and any pair of two photons (again 10 combinations per event). 
The four distributions for the data are simultaneously fit with the weighed sum 
of the same MC distributions for each background sample and for the signal. 
The weights of the rad and kk samples are the free parameters. Wrad — 0.45 
and Wkk — 1-3 are obtained, from which the numbers of background events 
Brad — 307 and Bkk — 264 are estimated. 8 additional background events from 
the all sample have also to be taken into account. The fit has been repeated 
separately on each control distribution and the spread obtained in the esti- 
mated number of events is taken as systematic uncertainty. The total number 
of background events is 579 ± 27, where the uncertainty is the quadratic sum of 
the statistical and the systematic uncertainties. This background accounts for 
about 14% of the selected events. 

FigO shows the rji:^ invariant mass distribution. In the same figure, the dis- 
tribution of the polar angle of the recoil photon is shown, and is compared to 
the MC expected behaviour. Also in this case the distribution agrees with the 
1 -I- cos^ Orec dependence of the signal. 
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M (MeV) cose 

T17I ^ ' rec 

Figure 5: Left: t^tt" invariant mass distribution for the final data sample (points) compared to 
the estimated background (dark histogram). Right: polar angle of the recoil photon for data 
(points) and for MC expectations (histogram). Dark histogram represents the background. 



3.3. Branching ratio evaluation 

The branching ratio of the process ^ ~> rji:^^ is obtained from the formula: 

Nt ~ Rt 

SfN^/^Brifj ^ /) 

where Nf is the total number of selected events, Bf the estimated background, 
£/ is the average efficiency. is the number of produced (f) mesons evaluated 
from the number Njj^ of — s- r/7 with r] tt'^tt^'tt^ events. 

Er^-yBrict) ?77)Br(?7 -> ttOttOttO) ^ ' 

The Br{TT^ 77) is not included in eq.(IT|) and ^ since it has been already 
taken into account in the MC. The normalization sample has been selected by 
requiring no tracks in the DC and six or more prompt clusters in the EMC, in 
the same runs used for the signal selection. TV^-y = 4.2 x 10^ events have been 
found in the sample used for the analysis of the fully neutral decay chain, with 
efficiency e^^ = 81%, corresponding to N^'^ = (1.24 ± 0.03) x 10^. 
By using Br{r] 77) — (39.31 ± 0.20)%[i0j, the branching ratio is obtained: 

Br{<j) 777r"7) = (7.01 ± 0.10 ± 0.20) x 10"^ (3) 

The first uncertainty is due to statistics and to the background subtraction. 
Several sources of systematics have been taken into account (see Tabled]): pho- 
ton counting (dominated by the detection efficiency for low energy photons), 
the data-MC discrepancies in the evaluation of the selection efficiency, and the 
normalization uncertainty. 

The data sample analyzed for the charged decay channel is slightly smaller than 
the other one, N^''~''°^ (1.15 ± 0.03) x 10^. By using Br{r] vt+tt-tt") = 
(22.73 ±0.28)%[y 

Br{<f) Tjn"-/) = (7.12 ± 0.13 ± 0.22) x 10"^ (4) 
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Table 2: Main sources of systematic uncertainty on the branching ratio !3l . 



Source 


uncert. (xlO ^) 


Photon counting 


0.08 


Selection efficiency 


0.12 


Br{ri 77) 


0.04 


Br{(j) — !■ r/7) 


0.13 


Brlri 7r°7r°7r°) 


0.05 



is obtained. The first uncertainty is the quadratic sum of the statistical uncer- 
tainty on N^+j^-T^o and of the uncertainty on the background; the second one is 
systematic, mainly due to the absolute normalization, and includes a 1% error 
due to the efficiency evaluation. 

The two branching ratios ^ and ^ are compatible with the old KLOE re- 
sults: (8.51 ±0.51 ±0.57) X 10-5 (?? ^ 77) and (7.96 ±0.60 ±0.40) x 10"^ (t] 
TT^TT^TT^) J>\. By combining the two results, taking into account the common 
normalization error 

Br((/) -> ry7r°7) = (7.06 ± 0.22) X 10-5 (5) 
is obtained, where the uncertainty is both statistic and systematic. 



4. Fit of the 7777" invariant mass distributions 

In order to extract the relevant parameters of the ao, a simultaneous fit, with 
the same set of free parameters, has been performed on the two rjir'^ invariant 
mass distributions, by minimizing the following x^- 



where is the number of bins of respectively the fully neutral and charged 
r^Tr" mass distribution; Ni is the content of the i-th bin and Bi is the number of 
background events to be subtracted from the i-th bin. The expected number of 
events, Ei^ can be written as 



bin j 



where m — Af,,^o, and — 4.26 MeV[lC|. e-^ is the efficiency matrix (also 
referred to as smearing matrix), representing the probability of a signal event 
with "true" mass in the j-th bin of the spectrum to be reconstructed in the i-th. 
bin. The efficiency matrices, evaluated by MC, are almost diagonal; the off- 
diagonal elements take into account resolution effects as well as wrong photon 



10 



pairings. The differential decay widtli dVth / dm lias been parametrized according 
to two different models. 

In the "Kaon Loop" (KL) modelfll] the is coupled to the scalar meson through 
a loop of charged kaons. The theoretical function can be written as: 

dm dm dm dm 

The scalar term dYscai j dm is described in some details in Appendix lAl dX^^ect /dm, 
is dominated by ^ — > pir'^ with p —^J2J and is described in the framework of the 
Vector Dominance Models fVDM')[12|. Last term is the interference between 
the scalar and the vector amplitudes. 

The free fit parameters are: the Qq mass, the couplings gaoK+K-i 9aovT^°i the 
branching ratio of the vector contribution, the relative phase 5 between scalar 
and vector amplitudes, and, as a relative normalization between the two differ- 
ent final states, the ratio Rj^ = Br{r] jj)/Br{ri 7r+7r~7r°). 
An alternative parametrization of the amplitude of the decay cj) 7/7r°7 has 
been also used, following ref.[l3|. A point-like coupling of the scalar to the (j) 
meson is assumed, hence this model will be called "No Structure" (NS) in the 
following. The scalar meson is parametrized as a Breit-Wigner interfering with 
a polynomial scalar background and with a vector background (see Appendix 
[B)) . The free parameters in this case are the couplings QaoK+K-^ Saorju-Oj and 
g^ao-yi the ratio the branching ratio of the vector background, and two com- 
plex coefficients, bo and bi, of the scalar background. The ao mass is fixed to 
avoid fit instabilities, due the large number of free parameters, and due to the 
large cancellations that occur among the terms of eq.©. The chosen value of 
the ao mass is the result of the KL fit. 

The fit results are shown in FiglGl and the parameter values are listed in 



Table 3: Fit results for KL and NS models. 





KL 


NS 


Ma, (MeV) 


982.5 ± 1.6 ± 1.1 


982.5 (fixed) 


9aoK+K- (GeV) 


2.15 ± 0.06 ± 0.06 


2.01 ± 0.07 ± 0.28 


gaov^» (GeV) 
g<t,ao'r (GeV^i) 


2.82 ± 0.03 ± 0.04 


2.46 ± 0.08 ± 0.11 




1.83 ± 0.03 ± 0.08 


S (deg.) 


222 ± 13 ± 3 




B.r. of vector backg. (xlO^) 


0.92 ± 0.40 ± 0.15 


- 




1.70 ± 0.04 ± 0.03 


1.70 ± 0.03 ± 0.01 


|bo| 




14.9 ± 0.6 ± 0.5 


arg(bo) (deg.) 




38.3 ± 1.1 ± 0.6 


|bi| 




21.3 ± 1.4 ± 0.9 


arg{hi) (deg.) 
X^/ndf 




57.3 ± 1.4 ± 1.1 


157.1 / 136 


140.6 / 133 


10.4% 


30.9% 
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M^^ (MeV) M^^ (MeV) 



Figure 6: Fit results: points are data after background subtraction; histograms represent the 
fit functions for KL (solid) and NS (dashed) models. 



Table 4: Correlation coefficients among the relevant ao parameters. 



KL model 


NS model 




Mao 


gaoK+K- 






9aoK+K- 




9<pao1 


Mao 


1. 






9aoK+K- 


1. 






9aoK+K- 


0.931 


1. 




9aor]7r^ 


-0.565 


1. 




QaoTj-n^ 


0.584 


0.550 


1. 




-0.138 


0.657 


1. 



Table [31 Good probability is obtained for both models. 
The ratio i?^ is in good agreement with the PDG value 1.729 ± 0.028[ic 
firming that the two samples are consistent with each other. 
A vector background smaller than the VDM predictions, (3 5) x 10^^ 12, 15l |. 



is found in both fits, indicating that the rjir'^j process is largely dominated 
by (f>^ aoj. 

In the KL case, the ao mass is in agreement with the PDG value (985.1 ± 1.3) 
MeV[l3l. A ratio of the squared coupling constants Rao — 9aoK+K- /daorj-^" ~ 
0.58 ± 0.03 ± 0.03 can be derived. The 5,^007 is not a free parameter of this 
model, but can be obtained according to the formula: 

(7) 

= 1.58 ±0.10 ±0.16 GeV"^ 

The ao width obtained from eq.® is TaoiMao) ^ 105 MeV. 
In Table [3] the correlation coefficients among the ao parameters are shown. 
The couplings QaoK+K- and gaori-n" of the NS fit and therefore the ratio Ra„ = 
0.67 ± 0.06 ± 0.13 are in agreement with the KL values. In the NS case (?0ao7 
can be determined directly and is compatible with the value of eq.([7|. From 
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this fit a total decay width Tag {Mat, ) — 80 MeV can be evaluated according to 
eq.dlOl). 

The systematic uncertainties on the parameters account for: (i) sensitivity to 
the fixed parameters (the uq coupling to ri'ir^, daov'T^"^ ^^'^ 9ipK+K- ™ the KL 
model, Mao the NS model); (ii) normalization uncertainty; (in) data-MC 
discrepancy of the fraction of wrong photon pairings (12% from data and 14% 
from MC). 



5. Unfolding of the J77r° invariant mass distribution 

In order to allow a better comparison with other experimental results and 
with theoretical models, the invariant mass distribution should be corrected for 
resolution and smearing effects. Therefore an unfolding procedure has been ap- 
plied to the rjTv'^ invariant mass distributions by using the method described in 
ref.[l3|- This is an iterative procedure based on the Bayes theorem, which does 
not require the inversion of the smearing matrix. 

The unfolding has been performed separately on both invariant mass distribu- 
tions before the background subtraction. The smearing matrices are the same 
used in the fits described in Sect|4l 

An initial distribution has to be provided as starting point of the iterative pro- 
cedure; the unfolded distributions obtained starting from the output of the KL 
fit or from a flat distribution in M,j^o differ by less than 3%. This difference 
has been taken into account in the uncertainty evaluation. 
The bin by bin average of the two unfolded distributions is used to calculate the 
differential branching ratio {1/T^){dr{(f> rjir^^) / dM^^o) reported in Table[5l 
The uncertainties are both from statistics (data and MC) and from systematics. 
The main contribution to the systematic error is the difference between the two 
unfolded distributions. The correlation of the contents of nearest neighbour bins 
of invariant mass is about 50%, for next-nearest neighbour bins is about 20%, 
and is negligible for bin distance greater than two. 

An additional uncertainty of 3% on the absolute scale has to be considered, 
according to eq.®. 

To check this procedure, the unfolded distribution has been fit to the KL model, 
without requiring any smearing matrix. The parameters values are in good 
agreement with those of Table [3] 

6. Conclusions 

A high statistics study of the process (j) ~^ ijtt^^ has been performed, by 
selecting the decay chains corresponding to 7] 77 and ry — > tt+tt^tt*'. 
Br{(l) r]Tr°-f) = (7.01 ± 0.10 ± 0.21) x 10^^ and (7.12 ± 0.13 ± 0.22) x 10"^ 
respectively have been measured. 

A simultaneous fit of the two invariant mass distributions has been performed, 
which shows that the two samples are consistent with each other. 
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Table 5: Differential branching ratio: m is the bin center, the errors are the total uncertainties, 
and the bin width is 6.35 MeV. 



lit 




lit 




CMpV"! 


V iV±C V J 


CMpV"! 

V iV±C V J 


\1\±XJ V J 






oo\j .00 


^.^0 ZIZ U.-LO 


uy ( .00 


n 1 Q _|_ 10 


9i,^f{ 71 

OOU. / 1 


9 qc _|_ n 1 4 

Zi.OO ZIZ U.-L4: 


704 94 


n 1 Q _|_ n 1 9 


kf\'^ OR 

OUO.UU 


9 97 _l_ n 1 

Zr.Z 1 ZIZ U.-LO 




n Q1 _|_ n 1 

U.Ol ZIZ U.IO 


RRQ 41 


9 qc _|_ n 1 q 

.^.00 ZIZ U.IO 


71 Q4 


U.OU ZIZ u.uo 


R7^ 7fi 


AO 4- n ^a 

Zj.'^Zd ZIZ u.xu 


79q OQ 


n OQ _|_ n 1 1 

U.OO ZIZ U.-L-L 


1 9 


9 CQ _L n 1 

^.Ocf ZIZ U.-LU 


79Q fi^ 


\J.OO ZIZ \J. ± 1 


RSS 47 
000. 1 


9 on _L n 1 4 

^.OU ZIZ U.-Lr± 


7"^^ on 


n CI _L n 1 

\J.O ± ZIZ u . ±0 


RQA R9 


9 Q9 _l_ n 1 Q 

Zi . ijZ, ZIZ U . J- 


749 


n CO _|_ n 
u.oo ZIZ u.uo 


Q01 1 S 
yui . 10 


q 1 Q _|_ n 90 
O.lO ZIZ u.zu 


7AR 71 

1 to. i -L 


U.U 1 ZIZ XJ.KJ 1 


Q07 


q q7 _L n 1 7 

0.0 1 ZIZ U. -L 1 


7^^ Ofi 


U .0 ± ZIZ U .U ( 


Q1 RR 


q 40 _L f) 1 7 

0.^0 ZIZ U . -L 1 


761.41 


0.94 ± 0.10 


920.24 


3.67 ± 0.17 


767.76 


0.99 ± 0.11 


926.59 


3.94 ± 0.17 


774.12 


0.99 ± 0.08 


932.94 


4.29 ± 0.25 


780.47 


1.08 ± 0.09 


939.29 


4.63 ± 0.25 


786.82 


1.30 ± 0.10 


945.65 


4.89 ± 0.21 


793.18 


1.27 ± 0.13 


952.00 


5.20 ± 0.22 


799.53 


1.42 ± 0.28 


958.35 


5.40 ± 0.28 


805.88 


1.63 ± 0.28 


964.71 


5.44 ± 0.33 


812.24 


1.71 ± 0.14 


971.06 


5.35 ± 0.22 


818.59 


1.79 ± 0.16 


977.41 


4.94 ± 0.21 


824.94 


1.66 ± 0.18 


983.76 


4.02 ± 0.19 


831.29 


1.82 ± 0.15 


990.12 


2.80 ± 0.27 


837.65 


1.96 ± 0.12 


996.47 


1.51 ± 0.32 


844.00 


2.13 ± 0.13 







Both models used in the fits, the </>— scalar meson couphng through the kaon 
loop (KL model) and the direct coupling (NS model), are able to reproduce the 
experimental t/tt" mass distribution. 

From the fit results that 4> — * ??7r°7 decay is dominated by ^ ao (980)7, since 
the vector contribution is very small, Br(e+e~ — > VP — > 7777^7) < 10~^. 
The fit allows also the extraction of the ao(980) mass and its couplings to 7777°, 
K+K-, and to the (j) meson. The mass agrees at one standard deviation level 
with the PDG value. The two sets of couplings obtained from the fits agree with 
each other. Using these couplings, a total decay width of the ao(980) in the range 
80 105 MeV is estimated. The ratio Ra„ = i^+i^-l 9^ ^ 0.6 - 0.7 is 
obtained. A large g^ao-y tias been found (1.6 -j- 1.8 GeV~^) suggesting a sizeable 
strange quark content of the ao(980). 
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A. Appendix: main formulas of the KL model[lll 

The scalar term of eq.® has the form: 

dTscai ^ 2\g^ji+K-g{m)\^P^^o{Ml - m^) 
dm ~ 3(47r)2M3 

where 



2 

gaoK+K-9aor]7T° 



The detailed formulation of the KL function g(m) can be found in ref.[ll|. 

Daoim) is the inverse propagator of the oq: 

Daoim) = Ml -m^+ ^[i?^en,b(M, J - Habim)] 

ab 

The sum is extended over all the possible two particle decays of the oq: ab = ?/7r°, 
K+K~, K°K°, and r;'7r°. 
The flo width is: 

r / ^ _ T,abl^^abim) _ E 

ab 

9aoabPab[m) 

ao [m) — — — {o) 

m iDTrm 

where: 



Pabim) 



'/ _ (TOg -I- mb)2 \ / ^ (mg - mb)' 



The parameters of the scalar term that are determined by the fit are the ao 
mass and the couplings gagK+K- and gao-q-n"- The oq to rj'ir'^ coupling is fixed 
either to gaon'^ro = -V2 cosippga„K+K~ iqqm hypothesis) or to gao-n'-^" = 
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2 siiKfipg^gX+x- [qq hypothesis), where ipp is the pseudoscalar mixing angle 
(the value tpp = 39.7° has been used[l7|). Another fixed parameter is the 
coupling of the (j) to the K^K~ pair: 



9,K.K- = (M| - 4Mi)3/4 = 4-49 ± 0-07 



B. Appendix: main formulas of the NS modelyj 

The diff'erential decay width of the NS model is the following: 



dm 3 A/| 



M2^+zM,„r,„(m) 

(9) 



The resonance width is mass dependent according to ref . : 

Tao (w) = r^TrO {m) + Tk+K- (w) + ^KORo (w) 

where : r^^o(m) = -^^Pr^^o; 

r^^M = 'Xm^' - mK/mf for TO>2M;^; (10) 

^KRim) ^ "Xl"' - 1 for m < 2Mk 

(with KK = K+K-, K^K") 

The scalar background is parametrized with a polynomial with two complex 
coefficients, bo and bi. The vector background, Ajj^ct, takes into account all 
processes e+e^ V ^ V Pi with V P27 {V,V' — p, uj, (f> and Pi 2 ~ 
ttO). 
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